Expansion of the Gibbs potential for quantum many-body systems: General 
formalism with applications to the spin glass and the weakly non-ideal Bose gas 
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For general quantum systems the power expansion of the Gibbs potential and consequently the 
power expansion of the self energy is derived in terms of the interaction strength. Employing 
a generalization of the projector technique a compact representation of the general terms of the 
expansion results. The general aspects of the approach are discussed with special emphasis on 
the effects characteristic for quantum systems. The expansion is systematic and leads directly 
to contributions beyond mean-field of all thermodynamic quantities. These features are explicitly 
demonstrated and illustrated for two non-trivial systems, the infinite range quantum spin glass and 
the weakly interacting Bose gas. The Onsager terms of both systems are calculated, which represent 
the first beyond mean-field contributions. For the spin glass new TAP-like equations are presented 
and discussed in the paramagnetic region. The investigation of the Bose gas leads to a beyond 
mean-field thermodynamic description. At the Bose-Einstein condensation temperature complete 
agreement is found with the results presented recently by alternative techniques. 

PACS numbers: 05.30.-d , 75.10.Nr , 05.30.Jp 



I. INTRODUCTION 

The problem of understanding the static properties of 
systems possessing large or infinite numbers of particles 
pervades all of theoretical physics. Apart from very rare 
exceptions, approximation must be employed to find the 
characteristic features of such many-particle systems. It 
is the static mean-field approximation which is usually 
used to find first descriptions of these many-body sys- 
tems. Although this approximation leads in many cases 
to reasonable descriptions there are situation where be- 
yond mean-field approximations are needed. 

More than two decades ago, the present author devel- 
oped a method [1] to derive such beyond mean-field ap- 
proximations in a natural way. The approach is based on 
the power expansion of the Gibbs potential and was per- 
formed for the Sherrington-Kirkpatrick spin glass (SK) 
model [2]. The investigation confirmed the Thouless- 
Anderson-Palmer equations (TAP) that had been ob- 
tained previously [3]. Due to the presence of infinite 
range interactions, the power expansion truncates. It are 
just the terms up to the second order which contribute 
in the thermodynamic limit. 

In subsequent publications, the method of [1] was suc- 
cessfully applied to the infinite ranged classical vector 
spin glass by Vulovic [4], to various other spin glasses 
[5, 6] and to dynamical problems in the field of spin 
glasses by Biroli and Cugliandolo [7, 8] . Neural networks 
are closely related to spin glass [9-11]. Therefore it is 
natural that the Gibbs potential approach has been used 
for problems in this and related fields [12]. In particu- 
lar the first, generally accepted TAP-like approach to the 
Hopfield model, worked out by Nakanishi and Takayama 
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[13], is based on the method of [1]. 

The study of non spin glass-like problems with the 
Gibbs potential approach started with the contributions 
of Georges and Yedidia for the ferromagnetic Ising and 
spherical models [14] and for the Hubbard model [15]. In 
addition to these investigations, further work for non spin 
glass- like problems [16] exists that uses the method devel- 
oped in [1]. In general, infinite range interactions are not 
presumed and therefore a truncated expansion implies 
an approximation and higher order terms may become 
important. For the general Ising model, the third and 
forth order terms have been calculated first in [14] and 
later in [13]. Moreover, for some special Ising or spherical 
systems, leading terms beyond the forth order have been 
presented in these papers. 

Note that the majority of the systems to which the 
expansion of the Gibbs potential has been applied are 
classical systems. Exceptions are the investigations for 
the Hubbard model [15], the fermionic spin glass model 
[6] and the quantum version of the spherical p spin glass 
model [8]. Apart from the fact that these approaches 
yield interesting results, none of them is completely rep- 
resentative for a generalization of the Gibbs potential ex- 
pansion to quantum systems. It is just the Lagrange pa- 
rameter conjugate to the order parameter and the chem- 
ical potential for which the Legendre transformation is 
performed in [15] for the Hubbard Hamiltonian. In gen- 
eral, the transformation is performed for a larger set of 
variables. Thus the approach is very special and there- 
fore non-generic for quantum systems. All operators - 
the spin operators, the number operators and the inter- 
action Hamiltonian - of the fermionic spin glass model 
[6] commute. Again such a situation is not representa- 
tive for a quantum system. Finally the approach [8] uses 
very special quantum variables and it is not obvious how 
to generalize this work to obtain explicit results for usual 
quantum spins. 
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Hence, it is the aim of this work to present a com- 
pete quantum version of the power expansion for the 
Gibbs potential. To work out the characteristic effects 
for quantum systems the approach should be as general 
as possible. Thus we start in Sec. II with an arbitrary 
Hamiltonian and work out all terms of the expansion up 
to an arbitrary order. In Sec. Ill the formal results are 
discussed from a general point of view. 

The remaining part of this paper is reserved for two 
specific applications. In Sec. IV the infinite ranged quan- 
tum spin s = 1/2 glass is treated for non isotropic inter- 
actions. TAP-like equations result which are discussed 
for isotropic interactions in the paramagnetic regime and 
compared to corresponding equations for classical vec- 
tor spins. In Sec.V the expansion of the Gibbs potential 
is worked out for the weakly non-ideal Bose gas up the 
first beyond mean-field contribution. The resulting equa- 
tions of state are valid in the entire temperature regime. 
For the self energy a complete agreement with previous 
and alternative approaches is found at the Bose-Einstein 
condensation temperature. In Sec. VI we present some 
concluding remarks and compare to other approaches. 

Some introductory remarks should be added on the 
physical systems to which the general results arc applied. 
The quantum spin s = 1/2 glass with infinite range in- 
teractions is the natural generalization of the SK or of 
the classical vector spin glass model to include quantum 
effects. Despite being formulated over two decades ago 
[17], an understanding of this spin glass model has proven 
elusive [18]. 

Although the theory of the weakly interacting Bose 
gas has a long history, some old problems - like the ef- 
fect on the critical temperature due to such interactions - 
have been revived (compare [19] for an overview). In this 
context detailed beyond mean-field investigations for the 
self energy have recently been published [20]. As these 
results were derived with both Greens function meth- 
ods and within the frameworks of Ursell operators [21], 
this approach is an ideal reference system to compare the 
Gibbs potential approach with other work. 

II. GENERAL FORMALISM 

A. Basic concepts of the Gibbs potential expansion 

The many-particle system is described by the thermo- 
dynamic Hamiltonian (/? = 1/T and fcs = 1) 

/C„ = -pn = uiA' + aJC' (1) 

i 

where the single particle contributions and the many- 
particle interactions correspond to J2i ^a^^ ^'^id to 
spectivcly. It is assumed that the thermodynamic quan- 
tities and, in particular, all expectation values can be 
calculated for the free Hamiltonian /Ca=o and that the 
problem is to find at least approximations for the inter- 
acting system. 



This is a very common question and typical realizations 
are: interacting systems of identical particles, classical 
or quantum spin systems and combinations of these sys- 
tems. The present approach to these problems is quite 
general and does not need any further requirements or 
assumptions. 

The parameters = (f^)* are real and the operators 
= (A^)^ are Hermitian [22]. The variable a represents 
an (in general complex) expansion parameter and the ly^^ 
exhibit a a-dependence which will be specified below. 
The Hamiltonian of the original problem corresponds to 
the value a = I. This implies that the values i^^^^i are 
given and fixed by the specific physical problem under 
investigation. Clearly this also implies that we have to 
set a = 1 at the end of the calculation. 

The index i is a shorthand index for both different par- 
ticles and different operators acting in the same subspace 
of the individual particles. For a Hamiltonian in second 
quantization the represent products of creation and 
destruction operators and the index i may become a pair 
index. 

The object is to calculate a thermodynamic potential 
that determines the relevant thermal mean values 

(...)^ = Tr...n^ (2) 

where the density operator TZa and the partition function 
Za are given by 

n„ = — and by Z„ = Tre'^" , (3) 

respectively. The usual choice for such a potential is the 
free energy which is proportional to InZ^. This work, 
however, is focusing on the Gibbs potential which is 
related to the free energy by a Legendre transformation 
and defined as 

G^{A') = In Z„ - Y,<^' ^itli = (■^')« • (4) 

Strictly, it is the quantity —(iGa{A^) which represents the 
thermodynamic Gibbs potential. We use, however, the 
term for Ga(A*), keeping in mind this difference. For 
systems with variable number of particles 7?.^ and Z^. 
represent the grand canonical density operators and the 
grand canonical partition function , respectively. This 
implies the slide modification /Ca = —jSlTi — jiM) for 
these systems where fi and M represent the chemical po- 
tential and the number operator, respectively. 
The total differential of Ga is given by 

dG« = (/C%da-^i/^d^^ . (5) 

i 

Thus the natural variables are the expansion parameter 
a and the variables which are conjugate to the La- 
grange parameters i/^. Thus both the Gibbs potential 
and the Lagrange parameters i/^ are functions of a and 
A^ and the present approach exclusively uses these nat- 
ural variables as independent variables. 
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Next Ga(A*) is expanded around a = keeping the 
values A* fixed to their thermal values. Thus the rela- 
tion A' = {A')a = {A')o = {A'}i holds. These con- 
strains determine the functions yl^{A>) and remove the 
arbitrariness of the above. The Taylor expansion leads 
to 



with G(") = 



5"G, 



Tl=l 



n! 



a=0 



(6) 



The zeroth-order term Sq{A^) of the expansion (6) is 
the entropy of the noninteracting system. This is eas- 
ily checked using Eq.(3), the general definition 5*0 = 
— (ln7?.o)o and 



(7) 



Considering next the first-order contribution G^^^ 
Eq.(5) yields immediately 



G(i) = (/C')o 



(8) 



which represents nothing else as the usual mean-field en- 
ergy. Thus approximating the expa,nsion (6) by the first 
two terms Ga=i ~ So+{IC')o corresponds to the standard 
mean-field theory. 

Differentiation of Eq.(6) with respect to and using 
Eq.(5) yields for a = 



and for a = 1 



v{=vi+ with - /JS' = E 



n=l 



1 aG(") 
n\ dA* ■ 



(9) 



(10) 



According to Eq.(9) the values of the Lagrange param- 
eters Vq are determined by the derivative of the entropy 
of the non interacting system that is a function of the 
Ai. This relation is important. Indeed, Eq.(9) is used 
to eliminate the dummy variables Uq which enter in the 
expectation values of the non interacting system via /Co 
(compare Eq.(7)). 

Recall that the quantities vl represent the given pa- 
rameter values of the original Hamiltonian (1) and note 
that both terms i^q and are functions of the A^ . There- 
fore the first Eq. of (10) represents a thermic equations 
of states from which the expectation values A^ can be 
determined for given values of the ul . 

The second equation of (10) represents the general def- 
inition of the self- energy or for magnetic systems the 
general definition of the internal field . Note that this fact 
is well known in the theory of Green's function where the 
Gibbs potential is often named effective potential [25] . 
Thus the expansion of the Gibbs potential implies a sys- 
tematic expansion of the self-energy that describes the 
effects of the interaction. 



For specific examples these results are illustrated below 
and the reader is referred to Sec. IV A and to Sec.V A for 
magnetic systems and for a Bose gas, respectively. 

In this context it is remarked that a further differentia- 
tion of Eq.(lO) with respect to Aj results in an expansion 
of the inverse susceptibility matrix. This quantity is of 
high importance as it governs the convergence of the ex- 
pansion and consequently the stability of the considered 
system - compare [1, 26] for the SK model, Eq.(83) for 
the quantum spin glass and Eq.(108) for the Bose gas. 

Rewriting the definition of the coefficients G^"^ as 

i 

it is obvious that the G^"' can be interpreted as a gen- 
eralization of the cumulant expectation values [23, 24]. 
It is just the additional a dependence of the Lagrange 
parameters i/^ that causes the difference to usual cumu- 
lants. 



B. Onsager term 

In this subsection the explicit expression for second- 
order term G^^^ = daaGa^o is investigated. This term is 
called the Onsager term and represents the lowest beyond 
mean-field contribution to the expansion. 

According to Eq.(5) daGa = (/C')c« holds which im- 
plies daaGa = Tr /C'(9tj7?.a • Focusing therefore on the 
calculation of daTZa the definitions 

Ea{U}= [ U{X)dX , U{\) = e^^^^Ue-^^" (12) 
Jo 



and 



(W|V)a = (W^Ea{V})c 



(13) 



are introduced. The definition (13) represents the well 
known Mori product [27] of the operators U and V. This 
product is a scalar product in the Liouville space, has 
many additional properties (compare Appendix 1) and is 
physically significant in particular for the linear response 
theory (see [28] for a general reference). 

The definition (12) permits us to represent the differ- 
ential rule for exponential operators as 

de'^° =E„{d/Ca}e'^° ; dJC^ = K.' da + ^Advl. 

(14) 

Due to the possibility of cyclic permutations within the 

trace Eq.(14) yields dZa = Za{AJCa)a ■ This leads to the 
total differential of the density operator 

d7^„ = E„{X;'}7^a da + ^¥.„{A}nc dvl (15) 



where U is defined by 

U=U- {U)oc 



(16) 
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The constrains = const, imply TrA''dTZa = 0. This 

leads to 

d.ul = -^r^iA^\IC% (17) 

j 

using the relation (A. 5). The matrix FJ^ is the inverse of 
the susceptibility matrix Xa 

E^"^" ='^^^- ^^^'^^ x^ = {A'\A^)a . (18) 

k 

As the Mori product is a scalar product the matrix Xa is 
positive definite and thus the inverse matrix exists. 

For a compact notation it is convenient to introduce 
the projectors by and Q„ 

F„U = {ima + Y,A'T'diA^U)a ; Qa = l-Pa- 

ij 

(19) 

Pq, and Qoc are super-operators which linearly map oper- 
ators of the Hilbert-space onto other operators of the 
Hilbert-space. With the above definitions it is easy 
to show that the usual projector relations P^ = Pq , 
Q2 = and PaQa = are satisfied. The projector 
Pa projects onto the subspace that is spanned by the 
elements and by the unit operator 1. These bases el- 
ements are linearly independent but are in general not 
orthogonal. 

Let us introduce some definitions. In accordance with 

[28] wc use the term observation level for the set of op- 
erators spanning the subspace Pq together with the 
unit operator. The set of all and the set of all are 
called constrained and conjugate variables, respectively. 
Using Eqs.(15),(17) and (19) the compact result 

5„7^« = E„{Q„/C'}7^„ (20) 

is found which directly leads to 

daaGc, = (/C'|Q„/C')a (21) 

and to the final expression for the Onsager term 

G(2) = (/C'|Qo/C')o . (22) 

Both the Mori product and the projector Qo are related 
to the Hamiltonian /Co of the free system. Thus G^"^^ can 
explicitly be calculated. According to Eqs.(7) the conju- 
gate variables i/g enter in this expression that must again 
be eliminated with Eqs.(9) to obtain the ^j-dependence 
of the Onsager term G^^^ . 



C. Cumulants for n > 2 

From the above treatment for the Onsager term it is 

obvious that higher derivatives are needed to calculate 
the G^"^ for general values of n. For this propose it is 



useful to generalize the definitions of E^, of Pq, and of 
the Mori product. 

First a commutative product Bi * B2 * ■ ■ ■ * Bn of an 
arbitrary number of operators is introduced [29]. The 
operation Eq on such a *-product is defined as a mapping 
to an (Hilbert space) operator given by 

E„ {Bi * B2 * . . . * Bn} = (23) 

/ dAi.../ dXnT[Bi{Xi)B2{X2)...Bn{Xn)] ■ 

Jo Jo 

The Afe-dependencc of the Bk{Xk) is given by Eq.(12) 
and T represents the thermodynamic (or imaginary time) 
ordering operator. It orders the Bk{Xk) operators with 
increasing Xk from the left to the right. 

Next the definition of the Mori product is generalized 
by 

(V|Si * ... * Bn)a = (V^E«{ei * ... * Bn})a ■ (24) 

The bra must always be an ordinary operator. It is just 
the ket which can be an ordinary operator, a ^-product or 
even linear combinations of these objects. Assuming in 
Eq.(19) that U represents such an object, the generalized 
projectors P^ and are still defined by these equations. 
The characteristic projector relation PaPa = Pa remains 
valid which again is easy to proof. 

The above generalizations imply new properties of the 
modified quantities. Some of these properties arc listed 
in Appendix 1. Two key properties, the derivatives of E^ 
and Pq,, are calculated in Appendix 2 and in Appendix 3, 
respectively. There we find 

daE„{Bi* ...*Bn}Tla= (25) 
Eq{(Qq/C') * Bi * ... * B„}TZa + Ea{da Bi*...* Bn}Tlc, , 

where the inner derivative 9„ Bi * . . . * B„ has be calcu- 
lated by the usual chain rule. 

Introducing the shorthand notation for the ^-product 
B = Bi * ^2 * . . . * B„ Eq.(25) and the definition (24) 
leads to the derivative of the generalized Mori product 

5q(ZY|B)q = (ZY|(Qq/C') * B)q + (5qZY|B)q + {U\d^B)^ . 

(26) 

In Appendix 3 it is shown that the derivative of Pq is 
given by 

da PqB = Pq (Qq/C') * (QqB) + Pq ^qB . (27) 

Both differentiation rules (26) and (27) are essential for 
the following. To this point it is assumed that B repre- 
sents an arbitrary ^-product. An extension, however, to 
linear combinations of such products is obvious. Assum- 
ing the usual addition, multiplication and differentiation 
rules for these linear combinations, the above Eqs.(26) 
and (27) also hold for linear combinations. Thus B rep- 
resents in general linear combinations of ^-products in 
these equations. 
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With all these extensions we arc well-equipped to cal- 
culate the higher derivatives of TZa by repeated applica- 
tion of the above rules . With the notation 



J^i — QolC' and 
we find 



1 



(forn>2) (28) 



with 

F(4 
F(5 



(a„)"7e„ = E„{Q„F(")}7ea 

= IC' 

+ 120J^2 + 120J^i * 



(29) 

(30) 
(31) 

(32) 
(33) 



(34) 



for low n values. Powers of the * multiplication are de- 
noted by bold power exponents (compare Eqs.(31)). 
For general values n > 1 the F^""*"^^ are given by 



f("+i) 



(n + 1) 



, k^\k2\...kn\ 

{ki,...,kn) 



T\^*Tr ■■■^n-' (35) 



where the sum runs over all fcj = 0, 1, 2, 
strain 



i=0 



n-\-\. 



with the con- 



(36) 



The general expression (35) can be proofed by mathe- 
matical induction. 

With these results we obtain for the cumulants 



= (/C'|QoF("))o 



(37) 



Note that Eq.(28) and Eq.(35) permit a recursive deter- 
mination of all the f'"^ Consequently this also applies 
to all cumulants G^"^. Some examples are 



= (/C'|Qo(Qo/C')^)o 



(38) 



= (/C'|Qo{(Qo/C')^ - 3(Qo/C') * [Po(Qo/C')^]})o 
= (/C'|Qo{ - 2(Qo/C')^ + 3(Qo/C') * [Qo(Qo/C')']})o 

which show a nested structure of the projectors. The dif- 
ferent forms of the G^^^ result just using the Pq + Qo = 1- 
All the expressions are treatable as they have to be cal- 
culated with respect by the bare Hamiltonian /Co ■ Again 
the dummy variables vl^ enter and have to be eliminated 
by use of Eq.(9). Thus all terms of the expansions for 
the Gibbs potential (6) and for the self energy (10) can 
in principle be calculated. 

Even for the classical system the representation of the 
higher cumulants based on projectors seems to be adven- 
turous compared to the approach based on generalized 



Maxwell relations [14]. Indeed relations like Eq.(35) for 
general n values have not been published for the latter 
approach. 

In this context it should be added that Eq.(38) applied 
to Ising systems leads to the known third and forth order 
terms for these systems [13, 14]. (For an explicit check of 
this claim the simplification (39) of below can be used.) 



III. DISCUSSION 
A. General remarks 

The expansion of the Gibbs potential (6) and the ex- 
pansion of the self-energy (10), together with the expres- 
sion (8) for the coefficient G^^^ and the relations (37), (28) 
and (35), represent the most general result of this work. 

According to the formal derivation the perturbation 
KJ is completely arbitrary. Moreover no restrictions en- 
ter for the bar Hamiltonian /Co or for the operators of 
the observation level. Thus one particle problems as well 
as many-particles problems can be treated with the pre- 
sented results. In the most general case the observation 
level may contain even many-particles operators. The 
work of Biroli and Gugliandolo [8] for the quantum ver- 
sion of the spherical p spin glass model represents such 
an interesting approach. 

It is always possible to extend an observation level by 
adding arbitrary operators to the original set. With 
the presumption that all the corresponding Lagrange 
multipliers 9^ satisfy the condition = the physical 
problem is not modified. Indeed the Hamiltonian and 
consequently the exact Gibbs potential do not change at 
all by such an extension. The projectors Po, however, 
differ which leads to different terms of the expansions. 
Thus different observation levels lead to different expan- 
sions. This implies that the quality of approximations 
like truncations of the power series depend on the chosen 
observation level. 

The latter conclusion is clearly illustrated by the fol- 
lowing simple limiting case. Adding formally the interac- 
tion A = KJ with = to any observation level Eq.(28) 
and Eq.(37) yields J^i = and G^"^^) = q , respectively. 
This implies that the first two terms of the expansion 
give the exact result whereas any truncation of the orig- 
inal expansion represents an approximation. 

The present approach requires that the bare Hamil- 
tonian can be represented as linear combination of the 
elements of the observation level [30] . Nevertheless a re- 
duction of the observables of the observation level is pos- 
sible. This can simply be achieved by transforming just 
a part of the parameters of the bare Hamiltonian to con- 
jugate variables. Again any reduction leads to different 
expansions. 

As far as no approximations are performed all obser- 
vation levels are equivalent. For specific problems this 

freedom can be used to choose a special observation level 
which leads to good or fast converging approximations. 
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Certainly such a procedure requires in general some phys- 
ical intuition for the system under consideration. Similar 
situations show up in all the approaches that are based 
on projector methods. In analogy it is generally expected 
that symmetry breaking operators should be included in 
the observation level apart from the bare Hamiltonian. 
For two particle interactions the complete set of all one 
particle operators should be an appropriate choice for 
the observation level (compare Sec. HID). Considering 
systems of Bosons or Fermions in second quantization it 
should usually be sufficient to span the observation level 
by all the occupation number operators (compare Sec. 
V). 

B. Classical systems and the quasi-classical case 



and 

cn = ((^'r)o (43) 

we find immediately 



(lnZ)(i) 


= Ci 




(lnZ)(2) 


= C2 


-c? 


(lnZ)(3) 


= C3 


- 3c2Ci + 2cf 


(InZ)W 


= C4 


- 4c3Ci - 3cl + 12C2C\ - &c\ 



These results agree with the Ursell-Mayer expansion, a 
cumulant expansion of the free energy for classical or 
quasi-classical systems [23] and demonstrates the close 
relation of the present work to these former approaches. 



For classical systems the * product reduces to the or- 
dinary product. All quantities commute and the Mori 
product (24) reduces to 

(V|Bi * B2 . . .)„ ^ (VtBiB2 . . 



•)c 



(39) 



which usually simplifies the calculation. 

For special quantum systems the situation may appear 
that all operators Ai and the interaction K' form a set of 
commuting operators. It is exclusively the algebra of this 
set which enters in the G^"-* and thus the replacement 
(39) can bo used to calculate these quantities. In the 
following we call such situations the quasi- classical case. 

The approach for the fermionic spin glass model by 
Rehker and Oppermann [6] is such a quasi-classical case. 
Thus it is obvious that the results of this fermionic sys- 
tem are very similar to the results [1, 3, 26] for the SK 
spin glass model [2]. Note that for this system a quan- 
tum treatment becomes necessary if a transverse mag- 
netic field is added. The existing approach, however, 
does not contain such a treatment. 



C. High-temperature expansion 

The limiting case of an empty observation level is of 
special interest. No variable at all is Laplace transformed. 
This implies a vanishing /Co and a quasi-classical situa- 
tion where Eq.(39) can be employed. The thermody- 
namic potential (4) is nothing else as the free energy In Z 
in high-temperature approximation. The statistical op- 
erator simplifies to 



no = l/Zo with Zo = Trl 



(40) 



and the general projector reduces to Pq which projects 
any U 



= {U)o 

onto the 1-direction. Setting 

(lnZ)(") 



n=0 



(41) 
(42) 



D. Consequences for two particle interactions 

In the typical many-body problem the interaction is a 

two particle interaction. Focusing on this case the ther- 
modynamic Hamiltonians is represented as 



/C„ = ^ I/; • A* + I ^ A' • C'^ with C 



0, 



(45) 

The indices i and j number the individual particles. For 
fixed i the components A^''^ and f of the vector matri- 
ces and represent one particle operators and the 
corresponding Lagrange parameters, respectively. The 
elements of the square matrices C'-' = C^' describe the 
interaction. To include the general case it is assumed 
that the observation level is spanned by the set of all 
(linear independent) single particle operators A*. 

The factorization property of the expectation values 
with respect to /Co leads to simplifications. The mean- 
field contribution (8) reduces to 



(46) 



Using Eq.(19) we find 

= i ^ 1^ . C''A^ with A' =A' - {A') , (47) 

which leads to the Onsager term 

G(2) = i ^ (1* • C'^ A^ 1 A' ■ C'^ A^)o. (48) 



The latter result implies several interesting features. 
First of all the Onsagc;r tc;rm is a superposition of the 
correlation functions which do not factorize in the gen- 
eral quantum case {i ^ j) 



{A 



^J,m,|^*,m.^j,m^-^^ (49) 
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In contrast, these correlation functions factorize for clas- 
sical systems or for the quasi-classical case. The presence 
of these quantum fluctuations can be the origin of essen- 
tial differences between classical and quantum systems 
(for an example compare below). 

Thermal averaging gives finite contributions to the On- 
sager term (48) only for such terms where each A"^'™' 
has at least one partner A^'"^^ . Thus the sums are dou- 
ble sums and triple sums do not appear in Eq.(48). In 
contrast to this behavior, the second order term of ex- 
pansions of the free energy InZ leads to triple sums as 
this term is given by (1|/C' * /C')o, according to [23]. 

A similar behavior holds for the higher cumulants. In 
any order the number of terms of the free energy expan- 
sion extends the number of terms of the Gibbs potential 
approach where at maximum n-fold sums arise. This 
conclusion is some indication that the correlations are 
more efficiently treated by the Gibbs Potential expan- 
sion. Note that in particular these arguments apply for 
systems with long ranged interactions. 

In this context it is remarked that a diagrammatic in- 
terpretation of the expansion can be given completely 
analog to the classical systems [14, 31]. As pointed out 
by these authors, the weak point of this method is that 
the vertex weight and the combinatorial factors can not 
be calculated systematically and thus the Feynman rules 
are not known. Nevertheless, some conclusions are pos- 
sible from these diagrammatic approaches. In particular, 
and relevant for this work, it is found that all diagrams 
are connected. Note that this can also be concluded from 
the general cumulant theory [23]. 

The last two conclusions cause the well known fact 
that for infinite-ranged models the expansions truncate 
for both non random and random interactions in the ther- 
modynamic limit. 

Consider first the case that the matrix elements of C"-' 
scale as N^'^ and are non random. From the discussion 
above follows all terms G^") with n > 2 are sub-extensive 
and can be neglected for large N. Thus the expansions 
Eq.(6) and Eq.(lO) reduce to the usual mean-field expres- 
sions 

G{{Ay) = So{{Ay) + l^{A)'.C'^A)^ (50) 

= C^'-N-')- (51) 

Recall that the function is the entropy of the non- 
interacting system. Thus the Eqs.(50) and (51) com- 
pletely determine all thermodynamic properties for the 
given parameters ul and C"-' of the Hamiltonian (45). 

Finally we consider random, infinite ranged systems 
where the matrix elements of C*-' are independent ran- 
dom variables, or where the C*^ arc proportional to a ran- 
dom variables. For these cases the scaling C*-' ~ N~^^^ 
has to be used to get the right N dependence for the 
extensive quantities. All cumulants with n > 3 are sub- 
extensive and it is just the Onsager (48) that must be 



added to Eq.(50). For classical systems these conclu- 
sions are well known and for quantum systems they are 
in agreement with [8]. 

IV. QUANTUM SPIN GLASS (s = §) 

A. General non isotropic case 

A system of A'" quantum spins Sj (with s = ^ and 
?i = 1) is considered in the presence of external fields 
hi- The spins interact via an infinite ranged spin-spin 
interaction and are described by the Hamiltonian 

n = -'^hi - Si-^'^Jij Si-Tsj , (52) 

where the dot denotes the scalar product in the three 
dimensional real space. The bonds — Jji ( with Ja = 
0) are independent random variables with zero means and 
standard deviations JN~^/'^. We consider a general spin- 
spin interaction and F represents an arbitrary symmetric 
tensor with real eigenvalues 'y^ and fi = x,y,z [32]. The 
norm of T is denoted by 

7' = (7^)' + (7^)' + (7^)' = trr^ . (53) 

The complete set of all one particle operators with 
i = I, . . . N and /i = x,y, z is used as observation level. 
Setting rrii — {si)a = {si)i and rui = \mi\ the entropy 
of the non-interacting system as function of the nii is 
well known and given by 

So = - + ™0 h\ + ^i) + - ™») • 

(54) 

The operator 

/Co = ^ t'i • Sj (55) 

i 

governs the calculation of the expectation values and the 
Mori products. Simplifying the notation the variables 
Uq and vl of Sec. II A are denoted by Ui and by (ihi , 
respectively. Prom Eq.(9) we find 

Vi = 2artanh(2rni) and — = — - = Si . (56) 

Vi mi 

Note that these equations have to be used to eliminate 
the dummy variables Ui. Eq.(lO) leads to the thermic 
equation of states 

Si 2 artanh (2m,) = [i{h^ - S,) . (57) 

According to Sec. (HID) all cumulants G*^"' with n > 3 
can be neglected in the thermodynamic limit. Thus the 
internal fields — Sj are given by 
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To calculate the terms G^^^ and G^^' we note that with 
Eq.(52) 



IC' 



Jij Si - T S -j 



(59) 



holds. All expectation values factorize, the mean-field 
contribution (8) becomes 



(60) 



and from the definition (19) 



QqJC' = ^"^Jij Si-T Sj with Si = Si -rui {61) 
results. 

The calculation of the Onsager term is straightforward 

but needs some more effort. Again using the factorization 
property and Eq.(A.7) the term (22) takes the form 



0- 



with 



(62) 



(63) 



This Mori product is treated in Appendix 4. The result is 
split in longitudinal, transverse and mixed contributions 



LT 



V'L/Li I \rL/l I V" 



+ ^ij , 



(64) 



which are given by 
1 1 



X. 



TT 



1 TO," 



1 (rrii + m 



(65) 

{{V%i-Tl) (66) 

^} tr (eiX)r(ejX)r 



m,: 



rrii 



Ui + 



+ 



rrii — n^i 



Vi - Vi 



} 



{7^-(r^),,-(r% + r?,} 



(67) 



The antisymmetric tensor associated with the cross prod- 
uct is denoted by (cjX). The quantities and (r^)ji 
are components 



(68) 



of the tensors T and F^. To write the Eqs.(66) and (67) 
as short as possible we have not completely eliminated 
the dummy variables which, however, can easily be 
done with Eq.(56) and with 



1 



dr 



(69) 



Note that the quantities l/vl and mi/vi have a physi- 
cal meaning. They represent the longitudinal and the 
transverse susceptibilities of the bare system. 

Putting things together we find finally with the usual 

replacement Jf^ J'^/N 



1,3 



and 



phi = Ui — P Jij Truj 



2N 



EdXij 
. dim 



(71) 



For given external fields hi the magnetizations rm are 

determined by the equation of states (71). Provided that 
these solutions nii are explicitly known all other thermo- 
dynamic properties follow from the Gibbs potential (70). 
Recall that these equation are exact in the thermody- 
namic limit. They are equivalent to the TAP free energy 
and the TAP equations for the SK model[l, 3]. 

To the bc^st of the authors' knowledge such results for 
the quantum s = 1/2 spin system have not been pub- 
lished previously. From the analogy to the TAP equa- 
tions it is expected that these equations will shed some 
light on the spin glass problem in quantum systems. Cer- 
tainly the present results are just the basis for this pur- 
pose and additional work is needed in this direction. 

Some elementary aspects are presented in the follow- 
ing subsections skipping points which are interesting from 
the spin glass point of view. Such questions need addi- 
tional efforts which are far beyond the scope of this work. 

Before we go into details a general aspect of our results 
is pointed out. Note that the transverse contributions 
(67) to the Onsager term exhibit energy denominators. 
Such denominators are a characteristic feature of all ex- 
pansions for quantum systems. The denominators are 
absent in the longitudinal and the mixed contributions 
given by Eq.(65) and (66) , respectively. As all operators 
commute this is obvious for the longitudinal part. For 
the mixed contributions it is a consequence of the cyclic 
property of the trace operation. 



B. Comparison with classical models 

Let us next work out the differences to the classical 
spin glass. For this we consider a system described again 
by Eq.(52) where Sf represent classical vector spins in 
three dimensions of length Sf — 1/2. Such a treatment 
leads to modified Eqs.(70) and (71). Obviously both the 
entropy term and the f j have to be replaced by the clas- 
sical expressions for S'q' and for the vf{m,i). 

As quantum fluctuations are absent the Onsager term 
simplifies (compare Eqs.(49)). We obtain 



{{s^j • rsf ) (Sf . r s^j))^ = tr xf rxf r , (72) 
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where the local susceptibility tensor xf is given by 

{xfY' = {{SfY(Sff)^ . (73) 

Writing again X^^ as sum of longitudinal, transverse and 
mixed terms we find for the transverse term 

» 3 

whereas the other contributions to X?- are still given by 
Eqs.(65) and (66) with Vi{mi) replaced by vf{mi). 

Note that no explicit representation exists for vf{rni). 
Nevertheless this function is well defined via the inverse 
function 



mi^^ coth 1^1 ^ 



and can be represented by the power expansion 

144to:^ 19008TO? 



(75) 



vf = Urrii 



5 +^^+0(mI). (76) 



It should be remarked that for the isotropic case F = 1 
all the results for the classical spins are in agreement 

with the previous work [4, 35] whore in addition the ex- 
plicit expressions for the classical entropy function can 
be found. 



C. Paramagnetic phase for isotropic interactions 



Employing the standard relations U = —dG/d/3 and 
S = G + /3U gives the internal energy U and the entropy 
S of this state 



U : 



-N^I3J'^ and S = N \xi2 - N ^ (3^ j"^ , (81) 
2,2 64^^'^^ 



respectively. For high temperatures /? ^ these results 
are expected from the high-temperatures expansion ac- 
cording to Scc.III C. In the low temperature regime, how- 
ever, both the internal energy and the entropy are not 
acceptable as both quantities diverge to the negative in- 
finity. Consequently there must be a phase with solutions 
rrii ^ for low temperatures in the zero field case. 

For further analysis we focus on the singularities of 
the susceptibility matrix xtj^ = dm^/dhj. The inverse 
matrix is determined by 

(X".!)^^ = dhf/dm^ = . (82) 

Reintroducing the expansion parameter a we find from 
Eq.(80) 

/3{x^r~^ = hi 5^-^ with It^ = a6ij -h^-f +^ (83) 



and with the coefficients 

5 {a(3J f 



a = 4-(- ■ 



12 



; 6 = a/3J ; c = . (84) 



From random-matrix theory [34] the eigenvalue spec- 
trum of the matrix lij is well known. It is a superposition 
of a continuous part 



To keep the discussion as simple as possible we spe- 
cialize to the isotropic case and set F = 1. This implies 
7^ = 3. To explore the paramagnetic limit m, ^ all 
terms have to be expanded. We find 



Ui = Arm + y m? + O (m?) 



and 



Xij = ^(5m2 + 5m2 + 2mi • m^) + O (mf ) . (78) 
lb 12 



(77) 



These expansions lead to the Gibbs potential 

G = Nln2 + N-^j3'^J^ for = (79) 
64 

and to the equations of states 

3 

+ ^(5m, + 7V-i5^m,)+OK) .(80) 

3 

Obviously the paramagnetic state (all rrii = 0) is a 
solution of the Eqs.(80) for vanishing fields hi = 0. 



X{x) = a — bx with 
and one discrete eigenvalue 



Xo = a + c + 



2 < X < 2 



(85) 



(86) 



which may be isolated from X{x). For the special values 
of the expansion parameter a{x) and ao the eigenvalues 
X{x) and Aq vanish , respectively. These values and their 
absolute values are calculated to 



a{x) = 2 



3a; ± i V 60 - 9^2 
50J 



ao = ±i 



^^^2 



(87) 



and to 



48 

|a(x)p = - (/3J)-2 |ao|' = 32 (/JJ)-^ , (88) 
5 



respectively. Realize that a{x) has a finite imaginary 
part for the possible a;- values ]a;] < 2 . Thus vanishing 
eigenvalues of the inverse of the susceptibility and sin- 
gularities for the susceptibility are only possible for the 
complex a values. In the complex plane the singularities 
a{x) are located on two sectors of the circle with the ra- 
dius |a(.T)|. As no intersections of these sectors with real 
axis exist there are no singularities for real values of a. 



10 



These singularities are of high importance, as aheady 
pointed out in [1], to which the reader is referred for 
more details . The tcrm-by-term treatment of the ap- 
proach of Sec. (II) can only be justified in the region of 
the complex a-plane in which the power expansion is 
convergent. The convergence criterium for a Taylor ex- 
pansion is given by \a\ < p where p is the radius of con- 
vergence. The distance from the origin (a = 0) to the 
nearest singular point determines this radius p . Thus, in 
the present case, p = \ce{x)\ holds [33] and the expansion 
for the paramagnetic solution can only be justified for 
temperatures T above the critical temperature 



(89) 



Tc = V5/3 J/4 . 



Below Tc the expansion and consequently the paramag- 
netic solution breaks down and at least one new phase 
must exist. 

For the SK model and the classical vector spin glass 
(compare below) the singularities of the susceptibility in 
the complex a plane show a similar behavior. For these 
systems the singularities are located everywhere on a cir- 
cle in the complex a plane. The circle intersects the pos- 
itive and the negative real axis, which implies two real 
values of a. From this point of view both the quantum 
spin glass and the classical models are quite similar. An 
important difference, however, results. For the classical 
models a physical staggered susceptibility diverges at Tc- 
Such a behavior is not found for the quantum spin glass. 

To complete the analysis of the paramagnetic phase we 
calculate the local susceptibility 



Xlo 



g{x) dx 



pJx + 5/12(/3J)2 



(90) 



of the quantum spin glass. Employing t he Wign er semi- 
circle law [34] for the density g{x) = ^(4 — x'^)/(2Tr) of 
the eigenvalues and introducing the reduced temperature 



(91) 



t = T/Tc = 4^3/5 T/J 



we find by an integration 



Xloc 



bif^ + 1)+ 4<2 - (t2 + 1) V25i4 - 10^2 + 25 

6J212 



(92) 

which holds for f > 1. The function xioc(i) exhibits a 
maximum at the critical temperature t = 1 and decreases 
with increasing temperature from the value Jxioc(l) = 
{i (7 - 2710) }i/2 = 0.474498 . We find xioc ^ /3/4 for 
/3 ^ 0, in agreement with the direct high- temperature 
expansion. 

The Eqs.(80) and (92) give a complete thermodynamic 
description of the paramagnetic phase above the critical 
temperature given by Eq.(89). Again the author is not 
aware of any work which has claimed these results before. 

It is of some interest to compare with the classical vec- 
tor spin glass. The classical results are well known [35] 



but can easily be rederived. We find with Eq.(76) by 

expansion 

In contrast to Eq.(78), terms proportional to • mj 
are not found. This results from the differences between 
{X^^r^ and 

The remaining calculation is completely analog to the 
quantum case and leads for the Gibbs potential to 



G"=' = const. + N{PJ)^ /192 



(93) 



The matrix I^j which governs the singularities is calcu- 
lated to 

(/(;.)^' = (12 + /32a2 J2/12) % - paJ.j . (94) 

The singular a- values are located on a circle with radius 
|a| = 12(/3J)~^ . This leads to a transition temperature 
of 



Tf = J/12 
and to a local susceptibility of 



(95) 



(96) 



for T > Tf . The results for all quantities differ. These 
differences are exclusively caused by the different equa- 
tion of states of the bare system given by Eq.(56) and 
by Eq.(76). The different expressions for Xf-^ have no 
influence. 

Obviously these conclusions are restricted to the para- 
magnetic phase at zero external fields. Incidentally we 
remark in this context that the 'quantum' term c/N of 
Eq.(83) is important for systems with an additional (in- 
finite range) ferromagnetic interaction Jq/N. Indeed in 
this case an additional contribution —Jq/N has to be 
added to Eq.(83) leading to a competition with the term 
c/N due to the different signs. 

The replica approaches [17, 18] claim for the critical 
temperature of the quantum spin glass Tc w J/(4\/3). 
These values are determined from the condition 1 = 
Jxioc(T'c) which is a consequence of the assumption [17] 
that near Tc the system has a continuous behavior, like a 
second order phase transition. Such a behavior, however, 
is not confirmed by the present work. Moreover Eq.(81) 
leads for a temperature of J/ (4\/3) to a negative entropy 
value of S" = —1.55685, which is impossible for a quan- 
tum system. Thus the replica theory values for Tc must 
be rejected. 

We conclude by noting that the discussion presented 
in this subjection can be extended to the general non- 
isotropic case and to the low temperature phase on bases 
of the presented results. For this propose all the tools 
[26, 36], developed for the understanding of the SK model 
on bases of the TAP equations, can be transferred to 
the quantum spin glass. Work in this direction for the 
quantum SK model in presence of a transversal magnetic 
field will be published separately [38]. 
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V. THE WEAKLY NON-IDEAL BOSE GAS 



A. The Gibbs potential in second order 



The model is described by the Hamiltonian 



2V 



q ' 



pqk 



where ek = k'^ /2m is the free-gas spectrum, the 6k, ^k 
bosonic operators, the A/k = ^k^k are occupation number 
operators and V is the volume of the gas. The interac- 
tion strength U is assumed to be momentum-independent 
which corresponds a delta function interaction in real 
space. 

We choose the set of all operators A/k as observation 
level. Thus the operators /Co, /Ci and JC' introduced in 
Sec. II are given by 



/Co=^i/kAAk , /Ci =^/3(/x-ek)Mc (97) 



and by 



pqk 



(98) 



Again we simplify the notation. The variables Vq and v\ 
of Sec. II A are denoted by z/k and by /3(/x — Ck) respec- 
tively. We work with the second quantization. Thus Ga 
represents the Legendre transformation of the logarithm 
of the grand-canonical partition function and /x is the 
chemical potential. 

Let us introduce the notation nk = (A/k)^. Then the 
entropy So of the noninteracting system as function of 
these variables rik is needed and explicitly given by 



Sq = ^{ (rik + 1) ln(nk + l)-ni, Inrik} 

k 

Employing Eq.(9), this leads to 

z/k = In rik - ln(nk + 1) 



(99) 



which has to be used to eliminate the dummy variables 
i^k- Note that an equivalent form of the latter equation 



nk = {exp(-t'k) - 1}" 



(101) 



is the Bosc function. 

The mean-field or Hartree Fock contributions are given 
by Eq.(8) and are immediately calculated to 




(102) 



where the relation (A. 36) is used. Similar to the spin 
glass system the treatment of the Onsager term is more 



complicated and therefore is presented in some detail in 
Appendix 5. Using the result Eq.(A.52) together with 
Eq.(99) and Eq.(102) we obtain the expansion of the 
Gibbs potential up to second order in U 

G(/J,nk) = ^{ (nk + l)ln(nk + l) -riklnnk} 

k 

f3U^ 

2^ npriq 

pq 

_^ ^ np+knq_k {1+ np+ riq) 



V 



I^^k ^P+k + ^q-k - i^P - i^q 



+ o{u^) , 



(103) 



where in principle all the i^k in the denominator of the 
second order term have to be replaced according to rela- 
tion (100). 

The general Eqs.(9) and (10) lead immediately to the 
equation of states 

/J(ek - M) = In ^ - - /^E^f ^ - ^E^^^ + O {U') 

(104) 

where E^"^' and E^^^ -t- E^^^ represent the first and the 
second order contributions to the self energy given by 

2U 



(1) Z,U 
^k = -y 2^"P ' 



(105) 



by 



^(2) ^ 2/3 [/^ ^ Ttp+qTlk-q - np(l -|- Tlp+g -|- Tlk-q) 



and by 



g(2) _2PIP 



^ (z/k + I^p - l^p+q - tlc-q)^ 



(106) 
(107) 



np(l 



tp+q 



)(l-hnk-q) (1 + np)np+qnk- 



(1 + nk) 



nk 



'(2) 

(100) Note that EJ^ results from the application of the chain 



rule for differentiation to (i/p+k 



^q-k 



Z^q)" 



The main part of the general results Eqs.(103)-(107) 
are well known in literature. Trivially this applies to the 
zero order contributions. All first order terms represent 
the usual Hartrce-Fock expressions [19, 20] for an inter- 
action that is momentum-independent. The second order 

(2) 

contribution EJ^ to the self energy looks like typical ex- 
pressions calculated by the Green's functions approaches 
[39]. 

Contributions like Ej^ ' that contain the square of 
(z/k + z/p — z^p+q ^ ^k-q) in the denominator are usu- 
ally not considered or discussed by the standard Green's 
function treatments. In this situation some arguments 
are presented for the relevance of this contribution in the 
next subsection. 
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Discussion 



The first argument for the importance of the con- 

— (2) 

tribution Sp is rather general. Note that this quan- 
tity enters in the susceptibility matrix that is defined as 
Xpq = —dnp/d{eq — n). Indeed the inverse matrix is 
given by 



Xpq 



1 d'^G 



(3 dn-pdn-p 



(108) 



Pripil 



2U 



as 



(2) 



One, 



This matrix governs the stability of the system and a 
fundamental property of this matrix is the symmetry re- 
lation = Xqp- Thus neglecting or modifying some 
terms of Xqp results in general to a violation of this re- 
lation with serious consequences. 

Various investigations exist for the weakly interacting 
Bose gas in literature. In particular the equation of state 
at the critical point of the Bose-Einstein transition was 
recently discussed in [20]. As argued in this work, the 
transition 



riq > 1 



(109) 



can be used in all equations. With this approximation 
we find in leading order 



(2) _ 2PU^ ^ np+qrxk-qnp [ ^ 



p q 



and 



"p+q "k- 



V(2)_ W 



Tlp-j-qT^k— q^p 



I "k _ »k Wk 

p q rip "p+q "k-q 



(110) 



(111) 



k = mode which becomes instable at the critical tem- 
perature. Such a behavior was assumed in the analysis 
[20] without any proof and the present approach confirms 
this assumption. 

Considering next the equation of states (112) we 
note that riq 1 implies the approximation (compare 
Eq.(lOl)) 



(114) 



where we have introduced in addition the 'dressed' en- 
ergy eq and the 'dressed' chemical potential jl as the 
q-dependent and the q- independent part of —Vq/ff , 
respectively. Setting the average value of the density 
n = J2p '^p equation of states (112) can be writ- 
ten in terms of the 'dressed' variables 



Ek — A* = ik —p, — 2Un 
2C/2 



(115) 



-I- 



7^ (^P+q - A)(ek-q - A)(ep - A) 



This is in complete agreement to the equation derived 
in [20], both within the frameworks of Green's function 
and of Ursell operators. Thus we conclude that the ex- 
pansion of the Gibbs potential is an alternative to other 
approaches for many-body systems of identical particles. 

Apart from this important conclusion for the present 
work we remark that detailed numerical investigations 
of the presented result are expected to be an interesting 
object of further research. Indeed, a complete numerical 
analysis of the Eqs.(103) and (104), which should include 
the entire temperature regime, may potentially give some 
new insight to the Bose-Einstein condensation for weakly 
interacting gases. This expecting is based on the impres- 
sive success of these equations for calculating the shift of 
the critical temperature of the transition [19, 20]. 



f'ands(^) 

energy are of the same order of magnitude, which demon- 
strates again the relevance of S^,^-* . 

Using riq ^ 1 we find as approximation for the equa- 
tion of states 



Both contributions T,\^ ' and SJ^ ' to the second order self 



Ck - M = 



2U 2/3[/2 

2^ "-P + -(72- 2^ "-p+q^k- 



(112) 

and as approximation for the inverse of the susceptibility 

Xkk' = ^ + — - -^72- 2^"q-Fk(2nq+k' + rik'-q) . 



I3nl ■ V 



(113) 

By a change of the summation index it is elementary to 
show the symmetry of Xkk' ' which is not surprising as 



both S^^^ and i:'^' are included in Eq.(113). 

All perturbations of Xkk' ^^'■^ order and can 

therefore be neglected [40]. This implies that it is the 



VI. CONCLUSIONS 

In this paper we have presented a systematic power 
expansion of the Gibbs potential for arbitrary many- 
particle systems including in particular all kinds of quan- 
tum systems. Developing and employing generalized pro- 
jector techniques we were able to present explicit formu- 
las which permit the calculation of the contributions up 
to an arbitrary order of the expansion for general sys- 
tems. After a detailed discussion of the general results 
the method is applied to two non-trivial systems, the 
quantum spin glass with infinite ranged interactions and 
the weakly interacting Bose gas. The; contributions up to 
the Onsager terms, which are the lowest beyond mean- 
field terms, have been worked out leading to new results 
or confirming recent results for these special systems. 

The present method has several advantages compared 
to other techniques. The method can be applied for all 
kinds of interacting systems including, in particular, sys- 
tems of identical particles, classical or quantum spin sys- 
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terns and combinations of these systems. No other tech- 
nique seems to have such a wide spectrum for appUca- 
tions. 

Compared to the efforts needed within the Green's 
functions or the Ursell operator approach the expansion 
of Green's is rather simple, direct and straightforward. 
In particular no partial summations are needed for the 
present approach to find the mean-field and the beyond 
mean-field contributions to the self energy. 

As a further advantage the present approach usually 
gives criteria directly for the convergence of the expan- 
sion. Within the framework of other techniques addi- 
tional investigations are usually needed to obtain this 
information. 

The application of the cavity method [9] is very com- 
mon for spin glasses and related problems. This method 
often allows convincing interpretations of the low order 
terms of formal expansions. The cavity approach was 
originally developed for classical systems. To treat in ad- 
dition quantum spin glasses work [41] has been presented 
which uses Trotter- Suzuki transformations and maps the 
quantum spin systems into classical spin models. Such 
treatments work only for special problems and can not 
be generalized to all quantum systems. Thus the existing 
extensions of the cavity method to quantum systems are 
restricted. The present Gibbs potential approach, how- 
ever, works for general quantum systems. 

Summing up we conclude that the power expansion 
approach may potentially represent a serious alternative, 
to the other, well-settled methods to treat the statics of 
many-particle systems. Certainly more applications must 
be worked out to confirm this possibility. 
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APPENDIX 
1. Some properties of the Mori scalar product 

Since the projector formalism is not very common, we 
list some elementary relations of this approach that are 
used in the present work. For further details and explicit 
proofs of these relations we refer to [28] . 

Apart from the general properties which are required 
for any scalar product the definition of the original Mori 
product (13) implies some additional properties. Let U 
and V be two arbitrary elements, then the relations 

(W|V)„ = (Vt|Z^t)^ = iy\uy^ = (z^t|Vt); (A.l) 
and the Kubo Identity 

{U\[lC^,V])o, = {[V,U^])a (A.2) 



can be proofed from the definition and the invariance 
of the trace to cyclic permutations. The property (A.2) 
implies the useful relation for special operators W. If 
[/Ca, W] = ujW with w 7^ is satisfied 

a;(V|>V)„ = ([W,V^])„ (A.3) 

results. 

The Eqs.(13) and (A.l) yield for the scalar product of 
the unit operator 1 with an observable 

(1|Z^)„ = (W)„ and {U\l)„ = {U^)„, (A.4) 

and for fluctuations of observables U = U — {U)a one 
finds 

(W|V)„ = {U\V)a = {li\V)a = {U\V)a - {U^) a{V) c ■ 

.... (^.-^^ 

The projectors and Qa sltc Hcrmitian in Liouville 
space and idempotent which implies with the definitions 
(19) the relations 

(Zi|P„V)a = {F„U\V)a = (PaZi|PaV)a (A.6) 

and 

(W|Q„V)„ = (QaW|V)„ = (QaW|Q«V)„ . (A.7) 

For later use it is noted that 

{A'^lQaKU = (A.8) 

according to Eq.(A.7) and QaA'' = 0. 

Let c be a complex number and let B = 81*82* ■■ ■*Bn 
be an *-product. Then the definition (24) immediately 
yields 

(W|B*cl)„ = c(Z^|B)„. (A.9) 
Again for later use we finally note that 

(l|(Qa/C') * (PaB))^ = (A.IO) 

holds. To prove this result we recall that PqB represents 
a linear combination of the unite operator and the A'', 
whereas Qa/C' does not contain such terms. Thus the 
thermal averaging eliminates all contributions. 

2. The derivative of Ea 

Let X be an ordered product of n operators Bk{Xk) 

A' = TBi(Ai)B2(A2)...B„(A„) (A.ll) 

where the A dependencies are given by Eq.(12). Assum- 
ing Ai < A2 < . . . A„ this product is already ordered and 

can be rewritten as 

(A.12) 
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Prom Eq.(14) one finds by elementary substitutions 



and obtains 



r'^'^(A)dA e^^+^^ 
Jx, da 



(A.13) 



+ B,{X,)J^ ^(A)dAB2(A2)...B„(A„) 



+ Bi(Ai)B2(A2)... / ^(A)dAB„(A„) 



A„-i 



5a 



provided that all the operators Bk are independent of a. 
The integral of the last term is rewritten with Eq.(12) as 



^ ^(A)dA= ^{x)dX-E^{dM- (A.14) 

Recalling the definition of ordering operator T the ex- 
pression for daX simplifies to 

daX = - XEc^id^K.^} + / dAT— -^(A)A'. (A.15) 

From Eqs.(14,17,19) 

i 

i 

= Qa/C' + C 

results, where c is a number. For such numbers the re- 
lations ]Ea{c} = c and TcX = cX hold. Replacing d^lCa 
by Qa/C' -I- c in Eq.(A.15) leads to 



3. The derivative of Pa 

Let B = Bi * B2 * • • • * be an ^-product. Then the 
definition (19) immediately yields 

a„p„B = a„(i|B)„ (A.18) 

+Y.^' r« ("^' 1^)" + E ^« ] • 

ij ij 

From Eq.(26) one finds 

(1|B)„ = (l|(Qa/C') * B)„ + (l|a„B)„ (A.19) 

and 

do. (I^|B)„ = (I^|((Q„/C') * B)„ + (I^|9„B)„ . (A.20) 

Again with Eq.(19) these two relations permit to rewrite 
the first two terms of Eq.(A.18) which leads to 

5„Pc«B = Pc,(Qa/C')*B + P„9c«B (A.21) 

Differentiation of Eq.(18) and employing again Eq.(26) 
leads to 

9a r^^' = - E ^» i^'lQalC' * A% r^' . (A.22) 

kl 

and with Eq.(19) to 

J2 [da ri] {A^ |B)„ = -J2K' {A'' I (Q„/C') * (P„B))^ 

j k 

+ J2 (11(Qa/C') * (l|B)„)^ . (A.23) 

k 

The Mori product in the second term of Eq.(A.23) can 
be written as (1|B)^ {A''\QalC')a according to Eq.(A.9). 
Due to Eq.(A.8) (^*=|Qa/C')a = holds and it is just the 
first term of Eq.(A.23) which remains. The multiplica- 
tion Eq.(A.23) with A^ and a summation yields 



daX = - A'E„{Q„/C'} + ^ dA T(Q„/C')(A) X (A.16) {A^\B)^ = -P„(Q„X;') * (P„B) , (A.24) 



as the terms proportional c cancel. 

Note that the result (A.16) does not change for any 
other order of the operators Bk in Eq.(A.ll). Therefore 
the above restriction on the A^ can be dropped and the 
Afe-intcgrations yields finally 

daEo.{Bi * ...*Bn} =Eo.{ Qa/C' *Bi*...*Bn} 

-E„{6i*...*6„}E„{Q„/C'} . (A.17) 

With Eq.(20) this result is equivalent to Eq.(25) in the 
text provided that all operators are independent of a. If 

the Bk depend on a it is obvious that the inner derivatives 
have to be added and one gets the full Eq.(25). 



where in addition the relation (A. 10) was used. 
Eq.(A.24) combined with Eq.(A.21) finally leads to 
Eq.(27) of the text. 



4. Calculation of Xij = (^s j ■ T s i \s i ■ F s j) ^ 

According to the definition (13) of the Mori product 
we have to calculate the quantity 



Xij = dx{{sj -Tsi} s,{X) -rijiX))^ . (A. 



25) 
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First the A-dcpcndence of Si{\) is considered. We set 
a = in the definition (12) and find from the rules of the 
spin s = I algebra 



(A.26) 



where the tensor r2j(A) describes a rotation of the imag- 
inary angle —i Af, about the axis 



(A.27) 



"r(A) 



-isinh(Ai'i)(eiX) — cosh(Ai/8)(ejX)^ 



and where 11^ and r2^(A) are the longitudinal and the 
transverse parts of fli{X) , respectively . The tensor 11^ 
is the projector onto the ej-direction and (eix) repre- 
sents the antisymmetric tensor associated with the cross 
product of two vectors. For further use the relations 



Hi 



and 



(A.28) 



nf(Ai + A2) = nJ{Xi)nf{X2) ; ni{X)mi = mi ; 
{fti{X)a}-b = a - ni{-X)b (A.29) 

are noted where a and b are arbitrary vectors. Using 
these relations and Eq.(A.26) the A-dependent part of 
Eq.(A.25) is rewritten as 

Si{X) -TsjiX) = {n,(-A)rni(A) Si} ■ Sj . (A.30) 

Let a and b be any two vectors (or two vector operators 
which commute with s j ) . From the well known identity 

{si ■ a){b-Si) = a • - ^(s,x)} b (A.31) 

it is elementary to prove the relation 

((s, •a)(b-Si))o =a- 0i(l/2)b (A.32) 
where we have introduced 



^~ +4cosh(i.V2) 



(A.33) 



and 



1 drrii 



i^-mi) . (A.34) 



Using the partial result(A.30), the relations (A.28) and 
(A.29) and applying the identity (A.26) twice we find 



= [ dAtr©; 
Jo 



(1/2- A) r ©J- (A- 1/2) r (A.35) 



where tr stands for the trace in the tree dimensional real 

vector space. The A dependence of ©i(A) is explicitly 
known and the integration finally leads to Eqs. (64-67). 



5. Calculation of Onsager term of the Boss gas 

First some elementary relations are deduced for later 
use. The definition rik = (A/k)o leads directly to 
drik/di/k = (A/'k)o - «k- From Eq.(lOl) we find 
dnk/dvk = nk(l + rik) and thus 

(A/'k)o = nk(l + 2nk). (A.36) 

Generalizing this procedure to higher order we find 

(^k)o = nk(l + 6 nk + 6 n^) (A.37) 

and 

(A/-4)o = nk{l + Unk + 36nl + 24 nl) . (A.38) 

Eq.(A.36) and the factorization property leads to 

{A\Afk')o = Skw {J^)o = ^kk' nk(l + rik) • (A.39) 

The Mc commute with ICq. Therefore {Afk\U)o = {■N'\M)o 
holds. This implies that the projector Pq defined by Eq. 
(19) simplifies to 



/O 



nk(H-nk) 



AAk 



(A.40) 



where U is any Hilbert space operator. 
Let us introduce a short hand notation by 



Spqk = hl^]J^\_^b^b^ 



(A.41) 



Then the interaction Hamiltonian (98) is rewritten as a 
sum of two contributions 



K 



with 



3;(i) = 2 E AApA/q + E [Ml N^] (A.43) 
3^(2) = Spqk . (A.44) 

pqk5^(0,q-p) 

With these definitions the Onsager term (22) can be ex- 
pressed as 

^^'^ = ^ {{ywi\ + {y\y2)^} , (a.45) 

where we have already used that (3^2)0 = and 
{Afj^y 2)0 = and consequently Po3^2 = holds according 
to Eq.(A.40). 

The two terms of Eq.(A.45) will separately be treated. 
As QoA/'pA/'q = AfpAfq results for p ^ q we find in conse- 
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quence that the commute with /Co 
(3^|Qo3^i)o = 2^(3^A/p7Vq)o + ^(mAA2)o 



p#q P 

p#q P 
8 ^ np(l + np)nq(l + riq) 
PT^q 

+4 5^n^(l + n^) (A.46) 



where the last step, the calculation of the commutator, 
is tedious but straightforward. 

With the definition (A.49) we find for k 7^ , q - p 



(3^|^pqk)o=4. 



^p+k rtq-k (1 + rip + Tig) 



^4 ripnq (1 + Tlp + k + »q-k) 



where in the last step Eqs.(A.36 - A. 38) are used. 

With Eq.(A.44) the second contribution to the Onsager 
term is written as 



{ym^ = E (3^l^pqk)o 

pqk#(0,q-p) 



(A.47) 



and we focus on the calculation of (3^|Bpqk)o for k ^ 
0, q — p . Note that 



[/Co , Bpqk] = E ''k [-A/ft , &p+k^q-k''p^c 



with 



Wpqk = t'p+k + fc'q-k - t'p - I'q 



■ Wpqk fipqk 

(A.48) 
(A.49) 



holds. Thus the relation (A. 3) can be employed that 
yields 

Wpqk (3^ I Bpqk )o = ( [fipqk, )o 

= 4([Bpqk,4qk])0 

= 4 { rip+k nq_k (1 + np + riq) 



where the Vp as functions of the rip are given by Eq.(lOO). 
Using this dependence we can calculate the limiting be- 
havior of (3^|Bpqk )o foi' the excluded values of k and we 
obtain 



lim (J^IBpqk )o = 4np(l + np)nq(l + riq) 
lim (>'|Bpqk)o = 4np(l + np)nq(l + riq) 

k^q-p 



for both cases. These findings imply that just one unre- 
stricted triple sum remains 



^(2) ^ W^U"^ r^p+k^q-k (1 + np + nq) 

pqk i'P+k + ^^q-k - i^p - I^q 



ripriq 



(1 + rip+k + '^q-k) }(A.50) and all the other contributions cancel out. 
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